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ABSTRACT 


Gravity anomalies have been recovered in the North Atlantic 
and the Indian Ocean regions. Comparisons of 62 2° x 2° mean 
free air gravity anomalies recovered in the North Atlantic area 
and 24 5° x 5 mean free air gravity anomalies in the Indian 
Ocean area with surface gravimetric measurements have shown 
agreement to *8 mgals for both solutions. Geoids derived from 
the altimeter solutions are consistent with altimetric sea sur- 
face height data to within the precision of the data, about 
±2 meters. 
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INTRODUCTION 

A spacecraft borne altimeter measures the height of the spacecraft above the 
instantaneous sea surface. In general, this surface devotes from the mean sea 
surface by no more than a few meters. The equipotential surface that corre- 
sponds to the mean sea surface (in the absence of dynamic effects such as ocean 
tides, currents, surges, etc.) is called the geoid. Since altimeter measurements 
are made relative to the instantaneous sea surface which closely follows the 
geoid, the altimeter measurements offer the best possibility for improving the 
accuracy of the marine geoid, and, by inversion, the marine gravity field. 

The GEOS-3 altimeter has demonstrated the capability to measure the fine 
structure of the mean sea surface. When this instrument is operating in the 
short pulse mode, a measurement is produced every 4 km footprint along the 
satellite subtrack. With this resolution capability it will be possible to describe 
sea surface topography to a detail of less than 1°, depending on the degree of 
data smoothing used and the spacing between subtracks. The data from the 
GEOS-3 altimeter constitutes an in-situ set of measurements of the sea surface, 
providing independent data over areas where surface measurements exist, and 
filling gaps in those regions of the sea surface where surface measurements are 
sparse or non-existent. 

In this investigation, emphasis is placed on the recovery of gravity anomalies 
from altimetry data in the N. Atlantic region bounded by latitude 20 J N to 40 “N 
and longitude 280 K to 200 E and the Indian Ocean area bounded by latitude 20’S 
to 50“S and longitude 90 E to 110 E. The recovered gravity anomalies are com- 
pared with gravity anomalies derived from surface gravity data only in the 
Atlantic area and with both surface gravity data and Apollo-Soyuz (Ref. 1' geo- 
dynamics experiment recovered values for gravity anomalies in the I i I . 

Ocean area. In addition, interccmparisons between recovered geoid, and altime- 
ter sea surface height and the 1® x 1 detailed geoid are studied for representa- 
tive altimeter passes in the Atlantic and Indian Ocean areas. 

The GEOS-3 altimetry data used in the solution for t hr North Atlantic region 
consisted of over 35,1)00 well distributed altimeter measurements obtained from 
portions of 80 altimeter passes which traversed the area. For the Indian Ocean 
area approximately 7200 well distributed GEOS-3 altimeter measurements were 
selected from over 72,000 measurements in 31 passes which traversed that 
area. 
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MEASUREMENT GEOMETRY AND MATHEMATICAL MODELS 

The geometry associated with the altimeter measurement is described in 
Figure 1. As can be seen, t' e altimeter is nominally the shortest distance be- 
tween the satellite and the sea surface, that is, the measurement is along the 
normal to the sea surface that passes through the satellite. 



The mathematical model for the altimeter measurement is given by the fol- 
lowing relationship 

h, = h - N' - h, - Ah' (1) 

where 


h S/C altitude above sea surface 

■ 


2 




■4 

h magnitude of h, the height of the spacecraft above the reference 
ellipsoid 

N' geoid height above reference ellipsoid 
h, deviation of sea surface from geoid 

h' systematic errors in altimeter measurement, e.g., refraction, an- 
tenna offset, timing, etc. 

The high rate (10/sec or 100/sec) binary format altimeter measurement data 
are processed in the GEODVN Program, (a general purpose o’* 1 . determination 
program which is used at the NASA Goddard Space Flight Center » where a tidal 
correction is applied and the data are converted to geoid heights and fit to a 
reference geopotential level surface. The tidal correction model in GEODVN 
includes corrections for both ocean tides (llendershott model) and solid earth 
tides (Ref. 2), but only the M2 component is considered. The GEODVN Program 
has been modified to permit the differential correction of orbital elements from 
altimeter data. This process can be used to converge the altimetric geoid 
heights to a profile through a reference geopotential level surface, in this case, 
the GEM-7 geoid (Ref. 3). The corrected orbit parameters serve merely as 
accommodation coefficients in this process and effect a removal of bias, trend, 
and long wavelength curvature from the relatively short arcs of altimeter data. 
While all six orbit parameters are adjusted, significant changes (greater than 
the standard deviation of the estimate) usually only occur in the orbit semi 
major axis, eccentricity and the true anomaly. Insignificant changes in the re- 
maining Keplerian parameters may also occur for the longer arcs. 

This process minimizes the square of deviation of the altimetric geoid heights 
from the GEM 7 geoid. 

To model the geoid surface, one must describe the gravitational potential. This 
is usually done »n terms of the potential of the reference ellipsoid, the normal 
potential, and the disturbing potential. The disturbing potential at the point 
P(0, \, r) is expressed in terms of spherical harmonics as follows: 


T 


GM v L i /a\ n n 


"10 L 


n a 2 


m*0 


[c n. cosm\ + S nm sin m\] 


( 2 ) 


where 


C 


n m 


G M 
S 

nm 


gravitational constant of the earth 
spherical harmonic expansion coefficients 
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P nrt (sill ;) associated Legendre functions 


(/', \) latitude and longitude at which disturbing potential 
is evaluated 

r geocentric radius from earth center of mass to evaluation 
[joint P 

a semi -major ;jtis of reference ellipsoid 

L is the limit of summation, and it is specified by the degree 
of harmonic expansion of the global geoid 

n summation index for degree terms of the spherical harmonic 
expansion of T 

m summation index for the order of terms in spherical har- 
monic expansion of T 

The geoid height at point P is expressed by 

N' N, ♦ AN, (3) 

where 

T 

N, - (Bruns' Formula Ref. 4), the global geoid model undulation 

y magnitude of gravity vector normal to reference ellipsoid 

AN, a local correction to the global geoid model values 

The global model contribution to the geoidal undulation at any point P(<£,\, r) 
on the geoid can be computed from geopotential coefficients derived from satel- 
lites by analysis of perturbations on the orbits induced by the Karth's gravity 
field. 

The Stokes' formula (Ref. 4, pp. 92-98) is another form of expressing the dis- 
turbing potential. This formula makes it possible to express the local details 
of the disturbing potential in terms of gravity data. That is, 

1 = ^1 f Ag S («/> ) d<7 (4) 

where 7 J*' 

R = a ( 1 ~ f) 1/3 

is the mean spherical radius of the earth, 


■ A 




f flattening of reference ellipsoid 

S(;) Stokes' function 

cr element of area 

£g surface gravity anomalies 
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(5) 


f rom Brum!' formula the undulation at any point P on the geoid car. b*. 

computed Irotn Stokes' formula (Equation '). That s* 


r ; r 

: — £(/ ) •» ’ 
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In terms of geographical coordinates, Stokes' function can be expressed as 
follows: 


N (sf . ) 


_r r ,n r 2 

4 • l: • J 


g ( X X* ) S( „ ) cos <t>' d <f>' d \* 
*V 2 


(7) 


where 


da cos d *>' d\' 


(/, \) latitude and longitude of the computation point 

(t‘ \') coordinates of the variable surface element * 

, spherical distance between the computation point and variable 
surface element 

* cos" 1 [sin<£ sin0* + cos <p cos cos(\ ~ \* )1 
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g( f>'^' ) free-air gravity anomaly at the variable point (/' >') 


y mean value of gravity over Earth 


To combine surface gravity data and geopotential information derived from 
gravity field perturbations acting on orbits of spacecraft, the Earth is divided 
into two areas (Ref. :t), a global area (A,) and a local area (A 2 ), surrounding the 
point P. Each gravity anomaly in each area is also partitioned into two parts 
represented by a g < and Ag 2 respectively, where the global model contribution, 
Ag s is given by 


A 


K. 



2>-»0 ».-< ^ ^ hr Ri ^ 

m*0 


(8) 


The Ag 2 value is defined as the local gravity anomaly. Uy partitioning the 
geoidal undulations into two corresponding components. Equation 7 can there- 
fore be rewritten as follows: 


N(4>. K) = N, + AN, 


(9) 


where 


N, 



Ag K‘) SO/-) cos <£' d<£' d>-' 


is the global model undulation term 


and 
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AN, 


R 

4r7> 



( <p\ \') S(^<) eos />'d*'dV 


ih the local undulation. Note that this partitioning also implies a partitioning in 
the disturbing potential. If we define 


and 

then 

or 

N 


AT 

AN, 

S(Ag t ) 


N(/.\) 


T ♦ AT 


N, 


RAg.*A\* 

4?’y 


b j A B, (4>i K\ )| S(V/,) cos 

i ■ l ‘ * 


( 10 ) 


where AN ( is the correction to the geoidal undulations of the global geoid as a 
function of the corrections of mean free air gravity anomalies. 

Kquation 10 is the form of the parameterization adapted for relating the altime- 
ter measurement residuals to geoidal parameters. That is by defining 


bN • Ah, AN, 

*8 * {Ag, <*J \* )} (11) 


A 


R A/* A\* 
Arty 


S(^ ) eos <p. 


Equation 10 is written in linear matrix form as follows: 



i) 


*<k ■ j > 


*R<i. 


>> 


( 12 ) 
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when* 


j - 1. 2, 3 m 

k »> j 


The gravity anomalies are estimated from altimetry data by using the weighted 
least squares estimation process. The solution vector i g ( ( . t) o( Kquation 12 
is obtained by solving Kquation 13 below. 



k 

-1 

■ k 

<* «>(,.») * 





.f*» 

(i •)> 

Lf-i J 


(13) 


where 

e 1 weighting matrix of altimeter measurement. 

W.'J a priori covariance matrix for anomaly blocks to be estimated. 

fcg 0 initial estimate to solution vector. 

The total gravity anomaly is computed by summing the local anomaly as ob- 
tained from (13) with the global contribution given by (9). That is: 

P, (<*>, \) ♦ 5 g (14) 

where 

* r {^r, (<fi. *•)} 


RESU LTS 

Area mean free air gravity anomaly solutions from altimetry data have been 
generated in the western North Atlantic and in the eastern Indian Ocean, west 
of Australia. Two degree by two degree equi-angular area means were esti- 
mated in the North Atlantic and 5 x y area means for the eas'.ern Indian Ocean 


8 



area. Kuch solution used the Goddard Karth Model-7 (GKM-7) geopotenti&l 
model for starting anomaly values. In « ach solution, the resulting geoid surface 
fits the ultimctrlc sea surface height observations closely, at or near the noise 
level associated with the altimek r data. This noise level, which varies with the 
changes in altimeter operating mode and sea surface phenomena, has been esti- 
mated by brown (Kef. 5) at between 1.8 and 2.1 meters for data in the North 
Atlantic. 

The quality of these solutions is difficult to measure since no accuiate stan- 
dard is yet in existence. Solution quality certainly varies from block to block, 
since the altimeter data density varies from block to block by one to two orders 
of magnitude. Aside from measures ‘>f self consistency (the fit of the solution 
geoid to the altimetric sea height data) the best test is to compare with other 
area mean free air gravity solutions. These other solutions are of two types: 
"ground truth" derived from point surface obsei vations of gravity and independ- 
ent solutions based on satellite to satellite tracking such as that obtained from 
the Apollo-Soyuz Geodynamics Kxperiment which used a different data type, and 
a different data processing technique than that described herein. 

Figure 2(a) shows the tabulation of 63 estimated 2° > 2' mean free air gravity 
anomalies in the Atlantic area and corresponding "ground truth" values based 
on 1 x 1 surface gravity measurements compiled by the Defense Mapping 
Agency, Verospace Center (DMAAC) The recovery region is bounded by lati- 
tudes 20 N to 40 N and longitudes 280 K to 300 K. For the solution, over 35,000 
altimeter observations from approximately 80 passes over the recovery region 
were processed while the altimeter was operating in the short pulse mode. 
Figure 2(b) snows the recovery region and the distribution of altimeter passes 
over that region. 

Dil.erences between "ground truth" and estimated anomalies for a 10 * 10° 
sublock solution are shown in Figure 3(a). Figure 3(b) shows the distribution 
of altimeter passes within this sublock area. 

Note in Figure 3(a), that large anomaly differences result for area blocks 1( 
and 22 which have the worst distribution of altimeter data and also art the 
largest in magnitude from the other anomaly differences. Upon their de’.ction, 
the RMS error is reduced from 12.4 mgals to 6.4 mgals. 

The accuracy ol the 2-2 gravity anomaly solution was checked by 

a) Computing a detailed gravimetric geoid using the 63 estimated 2° x 2 
mean free air gravity anomalies 
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GROUND TRACKS 

OF GEOS- 2 ALTIMETER PASSES USED 
IN SOLUTION 


(a) 2° X 2° ANOMALIES AND 'GROUND TRUTH 

(MILLIGALS) 


290 

E LONGITUDE 


3 4 5 6 7 8 4 290 " 234 56789 300 

E LONGITUDE 


RMS DEVIATION FROM GROUND TRUTH. 3.E MGA 


LEGEND: Q GROUND TRUTH'. DMAAC 

SPHERICAL DIST. 1^): 10° 

REFERENCE GRAVITY FIELD: GEM-7 


Figure 2. 2° x 2° Mean Free Air Gravity Anomaly Recovery Solution 
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Figure 3. Differences Between Estimated 2" x 2 ° Gravity Anomalies and Ground Truth 




b) Comparing the GEOS-3 altimeter geoid with the G8FC detailed geoid 
(Ref. 6) (Figure 4). 

The comparison of the GECXS-3 geoid with the detailed geoid for pass 1647 showed 
an unresolved offset of -2.1 meters and an RMS error of * 2 2 meters along this 
profile. 

From Figure 2(a), the RMS deviation from "ground truth" was found to be 
±9.6 mgal. The RMS error resulting from a comparison of the GKOS-3 altime- 
ter geoid with the detailed geoid is found to be *2.1 meters. Using the uncer- 
tainties cited above, a simplistic relation can be constructed which says that a 
RMS error in a 2 * 2 mean free air gravity anomaly of 4.0 rngals corresponds 
to uncertainty in geoidal undulation of • 1 meter. This value closely agrees with 
covariance error analysis results for which 1 meter in geoidal undulation 
equals *3 rngals for 2x2 mean free air gravity anomalies. 

GKOS-3 altimetry data taken over the East Indian Ocean have been analy/.ed to 
determine gravitatioi al features. Approximately 7200 well distributed altimeter 
measurements were selected from over 72,000 in 31 passes, each pass being 
individually fitted to he GEM -7 geoiu in the region bounded by latitudes 20 S to 
30 S and longitudes 90 K to 110 K. These selected measurements were used in 
a least squares estimation of 24 area mean free air gravity anomalies corre- 
sponding approximately to 3 • 3“ squares within the boundaries of the region 
(see Figure 3). The RMS error of commission for each of the estimated anom- 
alies was found to be on the order of *3 rngals, varying with the data distribu- 
tion within each individual block. Comparisons were made with 3 ■ 3 averages 
of 1 x 1 anomalies based on surface gravimetric measurements compiled by 
Prof. R. Rapp, Ohio State University as well as with independently estimated 
anomalies obtained from the analysis of satellite to satellite tracking (SST) 
Doppler data obtained from the Apollo-Soyuz Geodynamics Experiment of .July 
1973. In this comparison, most of the recovered anomalies were found to be 
within 0 rngals of Prof. Rapp's values with an RMS deviation of 9.9 rngals. For 
the comparison with the \pollo derived anomalies, agreement is found in all but 
two blocks to be within 8 rngals with an RMS of *7.9 rngals. Apollo derived 
anomalies also validate the GKOS-3 anomalies for those blocks where no 
"ground truth" values are available. 

A local marine geoid was computed from the 24 estimated gravity anomalies 
and this geoid was compared along two altimeter pass profiles w ith the GEOS-3 
altimeter sea surface 1 eight and the GSFC detailed geoid (see Figures 6-9). 

The two altimeter passes which are northeast to southwest (Figures 6 and 7) fit 
the altimeter geoid to within • 1.5 meters, whereas the southeast to northwest 
passes (Figures 8 and 9) fit the altimeter geoid to within 3.5 meters, in the 
former case, the altimeter data were much less noisy than in the latter case. 

In addition, SK to NW passes traverse over a larger portion of a major bathy- 
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Figure 4. Comparison of GEOS-3 Altimeter Geoid with Detailed Geoid and 
Altimeter Sea Surface Height Profile for Revolution 1647 
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Figure 5. 5 x 5° Gravity Anomalies Recovered from G EOS-3 Altimetry Data and Apollo/Soyuz SST Data 

in the East Indian Ocean Area 
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Figure 6. Comparison of GKOS-3 Altimeter Geoid with the Detailed Geoid and 
Sea Surface Height Profile for Revolution 1383 


metric feature, the Diamantina Trench, than the NE to SVV passes (see Fig- 
ure 10). 


CONCLUSION 

The results presented in this report should be categorized as preliminary, as 
more attention needs to be given to seeking the best way to treat orbit errors, 
tides, and deviations of the instantaneous sea surface from the geoid. It is also 
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(IN SECONDS) 

Figure 7. Comparison of GKOS-3 Altimeter Geoid with the Detailed Geoid and 
Sea Surface Height Profile for Revolution 1710 


to be noted that the "ground truth" solutions presented here are generally not 
of any better quality than the altimeter solutions. Formal statistics associated 
with these solutions are as large or larger than those of the altimeter solution. 

It can however, be concluded, that the feasibility of gravity anomaly recovery 
using altimeter data is proven. The altimeter data distribution available at 
present with the GFOS-ll satellite is sufficient for achieving 5" * 5 uiomaly 
block coverage over the oceans within the latitudes 65°S to 65 N. 
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Figure 9. Comparison of CKOS-3 Altimeter Geoid with the Detailed Geoid and 
Sea Surface Height Profile for Revolution 1718 
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Figure 10 



East Indian Ocean Area 
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